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Abstract 

A 1-parameter class of quadratic equations of state is confronted with the Hubble diagram 
of supernovae. The £t is found to be as good as the one using the standard ACDM model. 
However this quadratic equation of state precludes objects with redshifts higher than ~ 
1.7. Adding a fair amount of cold baryons to the model increases without spoiling the 
£t. 
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1 Introduction 


One attempt at easing the persisting tension between cosmic observations and theory is to 
admit an exotic matter component with an ad hoc equation of state expressing the pressure 
p of the component as a function of its energy density p\ 

P = fip)- (1) 

A popular example is the Chaplygin gas pQ and generalisations with /(p) = po + w p + a/p. 
Another functional class, quadratic equations of state, /(p) = po + w p + a p'^, has recently 
attracted attention: 

Using dynamical systems theory Ananda & Bruni |2] classify the many different be¬ 
haviours of the Robertson-Walker and Bianchi I universes resulting from quadratic equations 
of state. Linder & Scherrer [3] analyse asymptotic past and future evolution of universes 
with “barotropic fluids” i.e. matter with an equation of state ([^. 

Motivated by Bose-Einstein condensates as dark matter, Chavanis |1] has written a very 
complete series of papers on equations of state with /(p) = w p + a p^ in Robert son-Walker 
universes including analytical solutions of the Friedman equations and connections with in¬ 
flation. 

Bamba et ah [3] classify possible singularities of Robertson-Walker universes in presence 
of exotic matter, in particular with quadratic equations of state. Reddy, Adhav & Purandare 
[6] solve the Friedman equations analytically in Bianchi I universes with equations of state 
/(p) = —p + ap^. Singh & Bishi [8] consider the same setting in some f{R,T) modihed 
gravity theories. 

In the recent reference [7] Sharov confronts Friedman universes with curvature and general 
quadratic equations of state to the supernovae and baryon acoustic oscillation data. 

We would like to contribute to this discussion and concentrate on the particular case 
Po = 0 and w = —1. We justify our choice by Baloo’s mini-max principle: a maximum of 
pleasure with a minimum of effort. 

2 Friedman’s equations 

To start, let us write the Friedman equations without cosmological constant and without 
spatial curvature: 


3H^ = S7iGp, (2) 

2 H'+ 3 =-SttG p. (3) 

The prime denotes derivatives with respect to cosmic time t. We set the speed of light to one 
and have the scale factor a{t) carry dimensions of time. The Hubble parameter is as usual 
H := a'/a. We can trade the second Friedman equation (|^ for the continuity equation, 

p' = -3H{p + p). (4) 
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The coefficient a in the qnadratic equation of state has dimensions and it will be convenient 
to write it in the form, 


p =-p + =-p + I (5) 

where r is a characteristic time. Using this equation of state, the continuity equation Q 
integrates readily: 


P = 


1 

SttG 1 + 1 Hqt'^ ln(a/ao) ’ 


( 6 ) 


with integration constant Hq := H{tQ), the Hubble parameter today which is related to the 
density today and Oq := a(to) is the scale factor today, that without loss of generality can 
be set to oq = Is in flat universes. For vanishing r, we retrieve of course the cosmological 
constant A = p/iSnG). 

The pleasure continues and the first Friedman equation (|^ integrates as easily: 


a = ao exp 


{t - U)2/3 - [to - 

j2l-i 


(7) 


where is another integration constant. This scale factor agrees with the one obtained by 
Chavanis in appendix A of his second paper in reference |lj. 



Figure shows this scale factor together with the one of the ACDM universe, 

( cosh[A/3A (t - tbb)] “ ^ 

oacdm — 


ycosh[A/3A (to - tbb)] - 1 

with tbb the time when the big bang occurs. The fit is in anticipation of section 
Already now, the scale factor ([^ tells us a funny story: 


( 8 ) 
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• no violent, hot bang, just a little sneeze at t = tg- There, the scale factor remains posi¬ 
tive and the temperature hnite, only the Hubble parameter diverges together with the 
density and pressure of our exotic fluid. In the wording of reference [3] this singularity 
is of type III. 

• After the sneeze, the universe expands with deceleration until t± = ts + r where the 
deceleration parameter changes sign. Then the expansion accelerates for ever causing 
an event horizon. 

• The sneeze is a mild singularity, i.e. integrable and - if its hnite temperature is suffi¬ 
ciently low - transparent to light so that we can observe its past. There the universe 
contracts up to arbitrarily negative times, however with a particle horizon. 

• The scale factor is invariant under time reversal with respect to the time of sneeze tg. 

We are sufficiently intrigued by this universe to try and confront it with supernova data. 


3 Hubble diagram 

For the sneezing universe, the Hubble constant is 

2 (to - tgY‘1'^ 


H^= rr 


T 


5/3 


the redshift is given by 


, , ao (t- tgf/^ - (to - tgf/'^ 

z + 1 = — = exp--, 


and the apparent luminosity is 


(9) 


( 10 ) 


4iraJ V “o ’ 

with the absolute luminosity L and the dimensionless comoving geodesic distance 


f*to 


x{t) := / dt/a(t) 


3 

§ — exp 
ao 


[—x^/^exp -I- |-y/^erfx^/^] 


(t0-ts)/T 

{t-ts)/T 


( 11 ) 


( 12 ) 

(13) 


r 


Oq 


—R + I \/ti exp R^ erf R + exp 


t — t. 


r 


1/3 


exp — 


t -1, 


2/3 


— h \/n exp R^ erf 


t - tg 


l/3^ 


(14) 
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with the abbreviation 


1/3 


R:= 


(15) 


and the error function 


erfx : = 



exp Ay. 


(16) 


Note that the formula (14) is a priori valid only for emission times t G However due 

to the time reversal symmetry of the scale factor with respect to tg, equation (14) is valid for 
all t. In particular we have 


X{t) = 2x{ts) (17) 
Note that the redshift is not an invertible function of time, it has a maximum: 

I 



Figure 2: Hubble diagrams of sneezing and ACDM universes, colour codes as in £gure[^ 


'^max 


z{ts) = exp 


(^0 - 
^ 2/3 


(18) 


and takes negative values, blueshift, for t < 2ts — to- Its minimal value Zmin = —1 occurs at 
t = —oo corresponding to the particle horizon of finite comoving geodesic distance 


Xhorizon = x(-oo) = ^ R + exp R'^ {eiiR + 1}] . (19) 

Oo 

At this blueshift, Zmin = — 1, the apparent luminosity tends to infinity because the energy of 
the arriving photon is boosted during the long contraction phase. 

Figure 1^ shows the Hubble diagrams of the sneezing and the ACDM universes. 
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Because of the non-invertibility of the function z{t) the Hubble diagram is a genuine 
parametric plot containing two functions ^ 1,2 (^), 


n,2 


» = 


Qvrr^ (z + 1)^ -fp 2 (^) 

The function Ii(z) comes from emission times t G [tg, to) corresponding to z G (0, Zn 
Ii{z) = [ B? — \n{z + 1) dz 


( 20 ) 


( 21 ) 


= — R + \\fi exp W erf R 

+ (^ + 1) v/^R^^^MzTl)-|\/7fexpi?2erfv/)R^^^M^Tl) . (22) 

The function / 2 (^) comes from emission times t G (— 00 , tg] corresponding to 2 ; G (—1, ^^max], 

h{z) =2Ji(^max) — h{z) (23) 

= — R + exp erf R 

— ( 2 ; + 1) \/R? — ln( 2 ; + 1) + 1 1 /^ exp R? erf \/R? — ln(z + 1) , (24) 


where equation (23) is another manifestation of the time reversal symmetry via equation 

An example of a Hubble diagram with red- and blueshift, but without a cusp can be 
found in reference [9]. 


4 Sneezing universe versus supernovae 

To analyze the effect of a quadratic equation of state on the evolution of the universe, we use 
two data sets of type la supernovae. The hrst data set is from Union 2 sample [10] with 557 
supernovae up to a redshift of 1.4 and the second is from the Joint Light curve Analysis m 
with 740 supernovae up to a redshift of 1.3. The JLA analysis simultaneously hts cosmological 
parameters with light-curve time-stretching ag and colour at maximum brightness Pc- 

The following analyses use these 2 sets of supernovae independently and combined after 
statistical compatibility check. The combined sample contains 1007 independent supernovae, 
but since the light curve calibrations are different between the Union 2 and JLA analyses, 
two different normalization parameters are used. We use frequentist’s statistics [T2| 
based on minimization. The MINUIT package na is used to hnd the minimum of the 
and to compute errors by using the second x^ derivative. All our results are given after 
marginalization over the nuisance parameters {mg, ag and Pc)- 

The general x^ is expressed in terms of the full covariance matrix of supernovae magnitude 
stretch and colour including correlations and systematics provided by each collaboration. It 
reads 




( 25 ) 
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where AM is the vector of differences between the expected supernovae magnitude rUe and 
the reconstructed experimental magnitude at maximum of the light curve For the JLA 
sample the reconstructed magnitude reads: 


TTlr ^^peak (XgAl. 


( 26 ) 


where XI is related to the measured light curve time stretching, C the supernovae colour at 
maximum of brightness and rupeak the magnitude at maximum of the light curve fit. 

The expected magnitude is written as me{z) = rris — 2.5 log^Q £( 2 ;) where l{z) is given by 
the first branch of equation (20). Notice that the normalization parameters rus contain the 
unknown intrinsic luminosity of type la supernovae as well as the r parameter. The expected 
magnitude is then only a function of R and of the redshift (22). 

Table presents the results of the £t using Union 2, JLA and combined samples for a 
sneezing universe and the standard flat ACDM universe for comparison. The quality of the 
fits are identical for both models even though it is marginally better for the sneezing universe. 
This clearly indicates that the exotic fluid with the quadratic equation of state ([^ can replace 
both dark matter and dark energy at least for supernovae data. 



Sneezing Universe 

ACDM flat Universe 


R 

'^max 




Union 2 

0.99 ±0.03 

1.76 ±0.1 

530.78 

0.27 ±0.04 

530.73 

JLA 

0.99 ±0.03 

1.8 ±0.2 

701.07 

0.29 ±0.03 

701.29 

Combined 

1.00 ±0.03 

1.7±0.1 

1003.16 

0.30 ±0.02 

1003.18 


Table 1: Fit results (la errors) for Union 2, JLA and combined samples using flat ACDM 
and sneezing universes. 


Figure [^shows the data points and the black line is the best fit using only the first branch 
of the Hubble diagram (20). The red line indicates the second branch. We see that the red 


branch is not really populated justifying somewhat a posteriori that we have discarded it a 
priori. 

The parameter R for the combined sample is equal to 1 with an error of 3 % and using 


formulae (15) and (18) we have a maximum redshift of about 1.7, ^(max = e —1. Consequently 


the maximum temperature of the sneezing universe is 

Ts = (^max + 1) Tcmb = 2.72 ■ 2.72 K = 7.4 K, 
at the sneeze, some 9.1 Gigayears ago. 


(27) 
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Figure 3: The Hubble diagram for the three samples, Union 2, JLA and combined. The black 
line is the first, the red line is the second branch of the Hubble diagram computed from the 
sneezing universe. Only the black line was used in the £t. 


5 Adding baryons to the exotic fluid 

This £t encourages us to interpret the exotic fluid with its quadratic equation of state ([^ as 
ersatz of dark energy and dark matter together. We are therefore tempted to add ordinary 
matter, cold baryons with mass density pb and zero pressure, pb = 0. As we now have only 
one continuity equation. 


(P + f>b)' = -3ff (p + Pi+p). 


(28) 


















for two components, our system of differential equations is under-determined. To remain in 
business we cheat as is tradition by postulating two independent continuity equations for 
both components. These two continuity equations then integrate like charms yielding: 


flv : = 


Qh ■= 


8ttG p 


8 HI |i7oV2 1n(a/ao) + l/fixo’ 


SttG ph 

Tnf 


= fiw (- 


From the hrst Friedman equation ([^ we have for the initial conditions 

flxo + ^bO = 1- 


(29) 

(30) 


(31) 


Now remains the integration of the hrst Friedman equation (|^ with two components. The 
exotic huid alone generates the mild sneezing singularity, while the baryons alone generate 
the popular big bang singularity. 


a(t) 



- sneezing universe 

- sneezing universe + baryons 

---- ACDM 


Figure 4: Three scale factors 


With both components, the hrst Friedman equation has still separate variables and tells 
us that the cold sneeze kills the hot big bang: 


da 

a ^[( 9 / 4 )Ff 2^2 ln(a/ao) + 1/(1 - flxo)]”^ + fl6o(ao/a)^ 


Hodt. 


(32) 


Indeed going back into the past starting from to, the scale factor decreases until - at - the 
bracket [(9/4) Hq ln(a/ao) -|- 1/(1 — flxo)] vanishes. At this time, the scale factor has the 
value 


tts ■= a(ts) = ao exp 


4 

9 (1 - flw) HlH ’ 


(33) 
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Figure 5: Three Hubble diagrams, redshifts only, colour codes as in figure]^ 


and the universe sneezes and bounces. This is how it avoids the big bang. The maximum 
redshift. 


ao 1 4 

^max =-1 = exp —- .22 

as 9(1- Hbo) 


- 1 , 


(34) 


increases with respect to its value for vanishing baryon density. We performed the integration 


on the left-hand side of equation (32) numerically. 


Figure shows the scale factors for the ACDM universe, the sneezing universe and the 
sneezing universe with baryons added and hgure shows the three corresponding Hubble 
diagrams, redshifts only. Again the chosen parameters anticipate results of the £t to the 
Hubble diagram. 

Note that by the same mechanism, the sneeze continues to kill the big bang when we add 
radiation and spatial curvature. 


6 Sneezing universe with a pinch of baryons versus su¬ 
pernovae 

We use the same kind of analysis as in section 4 to study the sneezing universe now with a 
pinch of baryons added. The only difference is in the computing of the luminosity ^{z) from 


equation (20). As we have seen in the previous section, the evolution of exotic and baryon- 


fluid densities in terms of the scale factor are obtained from continuity equations and given by 
formulae (29) and (30). The hrst Friedman equation (|^ is solved with these two independent 
components using the Rugge-Kutta algorithm na with a step in time corresponding to an 
equivalent step in redshift well below the experimental redshift error (10“^) and by assuming 
again that we only see supernovae that exploded after the sneeze. The initial conditions are 
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a(0) = ao = 1 s, f2x(0) = and 1^6(0) = f2feo with the usual closure relation (31) for a flat 


universe. 

The hnal ht procedure is then only a function of time-stretching correction ag, colour 
correction (3^ flfeo and the new dimensionless parameter tHq describing the weight of the 
quadratic term in the equation of state ([^. Notice the simple relation between tHq and R: 
tHq = 2/3 ■ when flfeo = 0. 



^sneeze 

^baryons 

tHo 

R 

'^max 


Union 2 

1 n+'J- 

■‘-•U-0.15 

n H-U.lb 

U-UU o.66 

0.99 

1.77 

530.78 

JLA 


n iq+u-ib 

'^•■‘-^-0.14 


0.95 

3.2 

700.33 

Combine 

0.97;^;^^ 

n no+u.a 

U-U^_0,02 

n cv+u.ui 
U-O' -0.67 

1.00 

1.74 

1003.1 


Table 2: Fit results (la errors) for Union 2, JLA, and combined samples using the sneezing 
universe with some cold matter. 


Table shows our results of the £t for the three samples of supernovae. Despite the fact 
that we have one more free parameter, the quality of the £t is only marginally improved 
compared to a pure sneezing universe. However, the baryon density we hnd is statistically 
compatible with the generally admitted value of about 4%. The parameter R is very close 
to 1 as before and the maximum redshift is still quite small (of the order of 1.7) excepted 
for the JLA sample which gives Zmax = 3.2 because of the bigger amount of baryon density 
(14%). 

Figure 1^ shows the conhdence level contour on Dbo versus tHq. As expected, the degen¬ 
eracy between both fluids is important and explains the big errors in Table The best fit is 
indicated by the black dot. The numbers give the maximum redshifts. At a conhdence level 
better than 36%, both dark matter and dark energy can be replaced by the exotic huid with 
a pinch of baryons. 
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39% CL 
68% CL 
95% CL 


0.3 


0.2 


0.1 



Figure 6: Confidence levels for the sneezing universe with cold matter. The black dot repre¬ 
sents the minimum The numbers indicate the maximum redshifts Zma.^- 


7 Conclusions 

Just as dark energy with p = w p, the quadratic equation of state (|^ is a 1-parameter ex¬ 
tension of the cosmological constant. However the two ensuing universes differ dramatically: 
the second one has no big bang, except if r = 0. 

To fit the supernova data, dark energy needs the addition of cold dark matter and then 
retrieves the standard model. A nice feature of the sneezing universe is that it does not need 
cold dark matter to produce an excellent fit. Any amount of baryons, flfeO) may even be added 
without spoiling this fit. A further remarkable property of the fit in absence of baryons is 
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that with a precision of 3 % the characteristic time scale r in the quadratic equation of state 
(|^ coincides with the time span to — ts between the mild singularity, the ‘sneeze’, and today. 
We dually note that, thanks to its mild bounce without big crunch, the sneezing universe 
is an example of a “New Ekpyrotic Cosmology” [12], a possible alternative to inflationary 
models. 

However the sneezing universe has also dangerous properties: 

• It has no room for objects with redshift exceeding 

4 

^max = exp —- '22 “ 

9(1- Qbo) 

that is Zmax = 1.7 without baryons. 

• The temperature of the universe is always lower than its temperature at sneeze, 

(^max + 1) -Tcmb, that is 7.4 K without baryons. Therefore generically, it never ionized 
and was always transparent. 

• Without baryons the sneeze happened only some 9.1 Gigayears ago. 

• The sneezing universe predicts a cusp in the Hubble diagram at rather low redshift and 
a second branch coming from supernovae that exploded before the sneeze. This second 
branch is absent from the data. 

• This second branch also predicts supernovae with blueshift and high apparent luminos¬ 
ity. 

• Gibers’ paradox seems to make a come back in the sneezing universe. 

Certainly one can try and invoke new physics around the sneeze in order to avoid these 
dangerous properties, but we think that the sneezing universe deserves a cleaner dismissal. 

In order to confront the sneezing universe with Large Scale Structure data or lensing data, 
we need to know how the quadratic equation of state (|^ modihes the Schwarzschild solution. 
This question is currently under investigation. 

However we think that the cleanest dismissal comes from a confrontation with Cosmic 
Micro-wave Background data and we urge our colleagues practicing cosmic perturbation 
theory to take up this challenge. 
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